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ABSTRACT 

This paper consists to verify the validation of two multiple correlations between the 
Young ’s modulus (E) of a soil and its identification parameters. These correlations are used to 
calculate two types of theoretical moduli called Esikali and Elcpc ■ The exploitation of triaxial 
test results allows specifying 4 experimental moduli such as the initial modulus E t °, the 
equivalent modulus E e i q , the secant modulus at 1/2 and 1/3 of the breaking strain Eseci /2 and 
Eseci/s. Our approach compares the theoretical and experimental modulus through straight 
linear regressions. The results show that the Elcpc better represents E,° and E e i q while the 
Esikali is more consistent with Eseci/ 2 . 
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1. INTRODUCTION 

A joint study from the National Advanced School of Engineering (ENSP) and the National Civil 
Engineering Laboratory of Cameroon (LABOGENIE) has established multiple correlations between 
the Young’s modulus (E) and the parameters of identification of lateritic soils [13, 15, 22]. However, 
the Young’s moduli involved in the procedure has not been determined experimentally but deduced 
from empirical formulas. There are several different Young’s moduli, but empirical formulas proposed 
in the literature [20, 21] do not specify if the estimated E isinitial, equivalent or secant. 

This paper consists to verify the validation of those multiple correlations through experimental 
tests. For this purpose, the Young’s modulus is determined by adopting as soil behavior law 
hyperbolic model of Wong and Duncan. To obtain more specifications on the nature of the estimated 
Young’s moduli, the present study makes use of 26 samples of lateritic soils and the results are 
presented in order to show that the Elcpc better represents E t ° and E e i q while the E S ikali is more 
consistent with Eseci/ 2 . 

In this section, we introduce our work and section 2 deals with the materials and the experimental 
procedure of our approach. In this section we develop methodology that aims to verify the multiple 
correlations. The section III presents the different results and we conclude this study in section 4 with 
a short discussion. 

2. MATERIALS AND EXPERIMENTAL PROCEDURES 

To develop the methodology of our approach, we define the parameters of the multiple correlation 
model given in table 1 below: 


Table 1 Model parameters 


Parameter 

Description 

Mean of finding 

Drnax 

Maximum diameter of soil particles 

Sieve and hydrometer analysis [NF 
P94-056. NF P94-057] 

P 80 

Percent finer than 80 pm by weight 

D50 

Diameter of the particle at 50 percent finer 
on the grain size distribution curve 

G 

Percent by weight of particles coarser than 2 
mm and finer than 20 mm 

S 

Percentby weight of particles coarser than 20 
pm and finer than 2 mm 

L 

Percent by weight of particles coarser than 2 
pm and finer than 20 pm 

A 

Percent by weight of particles finer than 2 
pm 

PI 

Plasticity Index 

Atterberg limits tests [NF P 94-051. 
NFP 94-0521 

Elcpc 

Young modulus determined by the LCPC 
formula 

Elcpc = 5 CBR (E in MPa) [20] 

Esikali 

Young modulus determined by the SIKALI- 
MUNDI formula 

E<ukaii(E in bars ) 

(19CBR if P80 < 35%] 

- U 7 CBR if P80 > 35%) 
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The correlation models to be verified in this study are given in the table 2: 


Table 2 Correlation equations in [22] 


N° 

Power models 

1 

c — n0.94901 n/-0.32725 nO-07986 ^0.56787 rO.45942 /0. 13219 zi-0.29770 

£ LCPC ~ u max U 50 U b L A 

2 

Z 7 _ n0.79900 nj-0.36328 n 0.09940/-0.53986r0.39927 rO.10129 /1-0.18852 

£ SIKALI ~ ^mcix U 50 U b L A 


The procedure consists to compare the calculated and the experimental values of the Young 
modulus. In one side, by using the results of soil identification, the calculated value of E is determined 
with the formulas of the table 2. In the other side, the results of the triaxial test give the stress and 
strain values of the hyperbolic model of Duncan and Wong in order to find the experimental values of 
E. We examine four types of moduli, namely, the initial tangent modulus E t °, the equivalent elastic 
modulus E e i q , the secant moduli at half the breaking strain Eseci/2 and the secant modulus at one third 
of the breaking strain Eseci/ 3 . 


2.1. The hyperbolic model of DUNCAN and WONG 

In this model, the soil behavior law is based on an approximation stress-strain behavior of curves 
obtained by a drained compression triaxial test. The hyperbolic approach is the basis of Duncan and 
Wong model; the experimental curve [(g|- 03 ), Si] is approximated by the hyperbola of the equation: 

= (>) 

where a 3 , (a 1 — a 3 ) and s 1 are respectively the hydrostatic pressure, the axial stress and the 
vertical strain. The constants a and b depend on the triaxial test. 

The expression r ^ 1 = E^ is obtained by deriving the hyperbolic curve at the point with 

L ds-^ 

coordinate [s[, (eq — er 3 )i] . This is also the slope of the tangent or the increment of the Young's 
tangent modulus at this point. 

We have 

[for Ei = e[ 

[ E t = E‘ 

So E® correspond to = 0; E t ° is the initial tangent modulus. Then, we have 
E° = - (3) 

L a 7 

If we let (cq — a 3 ) uit — lim £i ^ OT [~“r] — ~ and the ultimate is denoted by the index ult , the 
hyperboliclawis translated as: 


riEi-oM 

dE-, 


2 - £ 

( a+be 1 ) 2 t 


( 2 ) 


Oi - er 3 ) = 


1 

“n + 


Ef (^1 ^3 ^ult 


We can also write: 


= Ef [l - 


a l~ a 3 
( Ci-a 3 ) u it . 


(4) 


(5) 
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(a) 


t (**“*») (o.-oju 



(b) 


Figure 1 Hyperbolic law: (a) represents E° and E[, (b) represents the law itself 


Wong and Duncan also propose to choose a point in the vicinity of the ultimate stress (asymptote) 
and set this point as breaking point; we then define a parameter Rf such that : 


(<Tl U3 ) ru p t 
( a l~ a -3)ult 


( 6 ) 


where the breaking is denoted by the index rupt and the ultimate is by the index ult. . 

Let C be the soil cohesion stress and <p'the angle of internal friction of soil under effective stress. 
The breaking point is also supposed to verify the Mohr-Coulomb criteria: 


We will have 


(oi - 03) 


rupt 


20 cos tp'+ 20-3 sin ipi 
1-sin (pr 


( 7 ) 


And 


(°1 °3 )nit — 


(<?i-G3) r upt _ 2C' cos <p' +203 sin <p' 
Rf Rf[l-sin(p'] 


( 8 ) 


E t = E" 


1 - 


R /(l~sin (p'')[(J 1 -rr 3 ) 
2 C' cos ^'+20-3 sin <p' 


( 9 ) 


It is often more convenient to effect a transformation of the hyperbolic pattern in linear diagram 
according to the following formulation: 


£ i _ 1 | £ i 

°i-o- 3 e? fri -^uit 


( 10 ) 



Figure 2 Linear representation of the hyperbolic law of Wong and Duncan [13] 
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The experience gained by Duncan and other mechanics of soil by performing simulations of stress- 
strain curves for different soils, showed that after the carrying of point of measurement onto the 
transformed linear diagram, these points are sometimes quite far from a perfect alignment. Thus a 
higher or lower concavity as far as the soil is dense or loose, is observed. 


2.2. Expression of experimental E based on Wong and Duncan model 


2.2.1. Expression of E® and E t 

According to Wong and Duncan [13], good setting can be obtained with the line: 

^ = /( £l ) (ID 

Passing through the two points corresponding respectively to 70% and 95% of the tensile strength 
(a x — a 3 ) rupt . In practice, these two points are carried into the transformed plane and is enough to 
define the right of transformation. 

The slope of this line is then given by: 


1 

(°1 — °3 )ult 


95% 

1 


(oi - 03) 


95% 

rupt 


,70% 


(ai - aj™ 


1 

f 95% _ £ 70% 


_ 1 /0.70£ 95% -0.95£[ 0% \ 

0.665( a 1 -cr 3 ') rupt \ £ 95% -£™ % ) 

where £^ 5% and £ 1 70% are respectively strain at 95% and 70% of breaking. 
We deduce 


1 

' p 70% -.95% ' 
t 1 t l 

1 

1 

1 

lO+j 

1 Iq 

„95% „70% 

)rupt 

( a l~ a 3)rupt. 


Then 


E? 2.66(o i - o 3 ) 


rupt 


,70% „95% 


.95% 


.70% 


(13) 


(14) 


From the expressions 


E t 


= E°t 



°1 a 3 
(°1 — a 3 )ult- 


and 


we obtain: 


( a l °3)ult 


(a i~(T 3 )rupt 2C 1 cos (p' + 2o 3 sin cp' 

Rf Rf (1-sin (pi) 



1 siny'frcrx-gg] 

2 C' cos cp' +2cr 3 sin cp' 


(15) 


(14) 


(17) 


2.2.2. The equivalent linear elastic modulus E etq 

The moduli are assumed reversible in the elastic area whether linear or non-linear. Consider 
respectively linear elastic law (Hooke's law) and the hyperbolic elastic law. £ e i denotes the linear 
elastic strain limit and £ eh , hyperbolic limits strainand (o 1 — o 3 ) el and (<r 1 — o 3 ) eh the corresponding 
stresses. 


As presented in figure 3a below for e\ < £ el , we have the linear elastic law: 
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Oi - 03 )e = E el 

In figure 3b for£{ < £ eh , we also have the hyperbolic law: 
(cti - a 3 y h = 


E° (tri-CT 3 )u/t 


(18) 


(19) 



Figure 3 Accent of the laws: (a) Linear elastic law and (b) hyperbolic law 
If one equates the hyperbolic law to Hooke's law we will obtain 

(0i-03 ')he= E el q £ i (20) 

With E e i q , being the equivalent linear elastic modulus and ( < o 1 — <T 3 ) l he hyperbolic stress 
equivalent to the linear stress. 

ACcq — <r 3 y, the error committed in equating the hyperbolic law to the linear law. 

A((7i - a 3 y = (ai - a 3 ) l h - (o-j - a 3 y he 

= (oi - 0 3 y h - E elq £ [ (21) 


The residual variation 


N 

S[A(0i - 0 3 ) ; ] = ^[(0i - 0 3 )h - E elq £[] 2 

i= 1 


N 

S[A (01 - 0 3 ) ; ] = ^ [[( 0 ! - 0 3 ))J 2 - 2E e i q £ l 1 {a 1 - a 3 y h + E 2 elq { sj) 2 ] (22) 

i=l 

To find the value of E eiq which minimize the error : 


We have 


dE e iq 


II- 

i= 1 


■24(01 - 0 3 )h + 2 F ei(? (4) 2 ] = 0 



(23) 


(23) 
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E elq — 


gfLiflCoi - o-iYh 

XIM) 2 


(24) 


With s u < e eh 

The expression of the incertainty on E e i q .This method is experimental, it occurs during testing 
difficulties affecting the measurement results: 

• The sudden closing of cracks at the beginning of change; 

• The friction of the piston ; 

• The imperfect sample size. 

The uncertainty in the determination of Young's modulus is then determined by differentiation of 
the natural logarithm applied to the expression of E elq : 


In E elq = In 


£f=iOi - °3) l h4 


= In 


- N N 

2 /oi - C? )h£i - ln[^(f{) 2 ] 

-i= 1 i = 1 


XIM) 

dE e i q _ d^jLiCgi ~ Ufi] _ d [Xdi(gj) 




E elq EfLiCOl - 

_ gLJdffil - + Oi - qgjhdfi] Xf= 1 2fid£l 

Sf = i(ai - Ef =1 ( f {) 2 

_ gyAffil - g 3 )ft£l + ffil - CTsUM] E"l 2£jAg{ 

A'eig SiliCof — °"3)h f l Ef = i(f{) 2 

Taking A[(cr 1 — (J 3 y h ] — 2A a and replacing in (25) , we have : 


(26) 


(27) 


A E e i q _ ZjLiJzAffehKffi-q^Aey _ 

E elq Eil 1 (o'i-o' 3 )J l £( Iili( £ i) 2 

Supposing As{ constant, (26) is rewritten as : 

A E eiq _ 2 A a g{ Az\ - cr 3 )£ 2 Ae) Iff =1 ej 

£ e ig - cr 3 )h4 + SfLiCo-i - o- 3 )he{ + Zf =1 ( £ {) 2 

Knowing that (crx — < 73 )^ < (04 — cr 3 ) e/l and£{ < s eh , we have 

N N 

_ ^ Oi - 03 ) eh ^ £ i ( 28 ) 

i= 1 i=l 

Considering 2f =1 (ffi ~ < 73 )^ ^ £ e /iEf=xOi - <7 3 )k an d ZjL i(^) 2 ^ £ e /i£f=i^ 
we will also have 


2 Act Z " =1 ej 


1 7 > 


2Acr 


Then 


- < 73)^1 (°1 _ °s)eh 

As[ Zf =1 (cr 1 - cr 3 )k Ae{ 


-7 — A > 


If=lOl “ <73)^1 
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and 

2As[ 27=i £ i 2Ae[ 

sr= i(fi) 2 £eh 

where 


(31) 


A E eiq > 3As[ 


2Ao 


J elq 


£ eh (°1 ~ a 3) 


eh 


We consider the expression of s leh et ( a x — <J 3 ) eh . The elasticity condition is given by the 
formula /[(ci — rr 3 )^J < 0 where f\{o l — <r 3 ) L h \ is the load function limiting the elastic range. 


/[(cr 1 - cr 3 )h] = Oi - ^3) ft - Oi + cr 3 )ft sin - 2C' COS( p' 
We will have at the limit of the elastic range: 

(ai - cr 3 ) eh - (d! + d 3 ) eh sin @ - 2C' cos <p' = 0 
2[C' cos @ + d 3eh sin@] 


(32) 


(°1 ~ ^3) eft — 

But we know that 


1 — sin@ 


(33) 


Hence 


(°1 — ^3) eft — 


£ eh 

1 | £ eh 

E? (o-i-o- 3 ) u ( t 


£ eh 


Oi - Os) 


eft 


70 

't 


^ (°~l~ cr 3)e/t 

(o-i-o-siuit- 


2 . 2 . 3 . Expression ofE secl/3 and E secl/2 

q * 

Generally, E sec 1/2 , secant modulus corresponds to deviator y by writing : 

g 50 _ f J_ 

Oi - cr 3 ) (d! - d 3 ) uit £s0 E° 

It therefore comes as: 

£ S eel/2 = E? (l - y) 

It is the same as E secl/3 — E° (l — y) 


(34) 


(35) 


(36) 


(37) 
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Figure4. Representation of Esecl/2 and Esecl/3 
Off (°1 — °3)rupt 9a — (@1 — °3 )ult) 


2.3. Experimental Determination of Parameters 

2.3.1. Presentation of Parameters to be Determined 

By modeling the law of behavior of lateric soils using the hyperbolic laws of Wong and Duncan, the 
determination of Young Modulus and their uncertainty necessitates the measurement or the knowledge 
of the following parameters for every test j : 

• The hydrostatic pressure er 3; - 

• The axial stress or deviator (aj — of) j 

• The vertical deformation £ 3 j 

• The soil cohesion effective stress C’ 

• The angle of internal friction of the soil under effective stress <p' 

• The deformation at 95% of breaking £^ 5% 

• The deformation at 70% of breaking £ 1 70% 

• The deviator at 95% of breaking 95%(<x 1 — o 3 ) rupt 

• The deviator at 70% of breaking 70%(n 1 — o 3 ) rupt 

• Coefficient Rf 

• The hydrostatic pressure in the hyperbolic yield o 3eh 

2.3.2. Basic Information 

The experimental determination of the parameters was carried out mainly on the basis of triaxial test 
results by LABOGENIE of 26 soil samples including 11 samples from Bini A Warak site and 15 
samples of LOM PANGAR site. The device used is of the "Stady" mark in Figure 5. 
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Figure 5 Photograph of a sample instrumented left and schematic view of the statistical measurement system 

(local) and dynamic (radial) 

The principle of triaxial test consists of a cylindrical test sample a certain number of loading cycle 
which are the resultant of lateral pressure or confining variable cr 3 pressure and axial stress u 1 equally 
variable. Displacement sensors allows for the measurement of the axial deformation (fq ) (Figure 6). 
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£ 




£ 


£ 


7 cm 


Figure 6 Representation of pressures applied to a test piece of the triaxial test 

In practical terms, we fixed er 3 and gradually increase o x until rupture of the specimen. This will 
permit us to draw for a given value of er 3 the curve (a 1 — cr 3 ) = /(fq). Tests in LABOGENIE were 
performed with three values of cr 3 (lOOkPa, 200kPa and 3 00k Pa). Thus, each soil sample is 
characterized by three stress-strain curves. Figure 7 shows the representive curves of the sample N°5 
fromBini A Warak site 
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Figure 7 Representatives curves (o j — a 3 ) = tfE}) of the sample N°5 from the Bini A Warak site 


2.3.3. Determination of parameters C ’ and <p' 

They are the parameters of elastic -plastic model of Mohr-Coulomb and these are determined at the end 
of the triaxial test. Thus, after drawing the stress-strain curve, we take for each level of confining 
stress, the maximum value of the stress deviator (effective stress). Then we place the point of abscissa 

t = ° 1 2 ai ‘ and s = in lame plan (s, t). Hence we get 3 points in the plane corresponding to each 

level of confinement. Using the linear regression method, we introduced a line called regression line, 
ensuring that the gap between the line and the various points is minimal. 

We know that in the elastic range, we have: 


°1 °3 °i + °3 . , , 

■sin cp + C cos cp 


(38) 


2 2 

which is a line ( y — ax + b) in the Lame plan 
Therefore by identification, we have: 

• The slope is a = sin cp’ 

• The intercept is : b = C' cos cp ' 

We have realized for the sample N°6 of Bini A Warak site the calculation of the parameters and 
the next figure below presents the results obtained: 


3.5 
3 

2.5 
2 

rsl 1.5 


y = 0,938493023x ♦ 0.04488 766, 
R* = 0.938 


CD 


1 

0.5 

O 


1.684 2-6435 

(cj’1 + cj' 3 )/2 (KPa) 


Figure 8 Regression line corresponding to the sample N°6 of the Bini AWarak site 
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Using the regression line, we get the parameters as it follows: 

sin cp' = 0.938493023 => <p' = 34.9 

0.4488766 

C' cos <p' = 0.4488766 => C' = = 13 kPa 

cos cp 

The procedure is applied to all 26 samples and the results are presented in the following paragraph. 


2.3.4. Determination of 

parameters 95%(<r 1 - (r 3 ) rupt , VOVoCa-f - <r 3 ) rMpt , £™ % , £^ 5% , Rf, and a 3eh 

The value of the deviator stress at breaking point is the maximum value of deviator. It is therefore 
trivial to calculate the values of 95% (a x — a 3 ) rupt and of 70%(a 1 — a 3 ) rupt . In regards to eI° % and 
e? 59 \ using linear interpolation with two points flanking points respectively the points[95%(a 1 — 
a 3 )ruptl £ i 95 %] ar| d [70%(a 1 — o 3 ) rupt ; £^ 0% ] one can readily determine the parameters eI° % and 
£ 95 % j|^ lcsc 4 parameters allow us to determine the value of the ultimate stress: 

1 1 / 0.70£^ 5% - 0.95£i° % 

Oi - cts )ult 0.665 (o'! - er 3 )™ P t \ fi 5% - el° % 

As well as 

1 1 f fl 70% £l 95% \ 

E° t ~ 2.66 (o'! - a 3 ) rupt [e™ - £ 70% ) 

The parameter Rf is the ratio between the breaking deviator and the ultimate deviator. To simplify 
the calculations, we used Excel software to get the results from the N°1 sample from Bini A Warak 
site [table 3]: 



Table 3 Parameters of the model of sample N°1 


Stress 

Strain 

Initial secant 

modulus 

Ultimate 

stress 

Coefficient 

03 

(kPa) 

(°1 a -i)rupt 

(kPa) 

(°i ^ 3 ) 95 % 

(kPa) 

(°i ^ 3 ) 70 % 

(kPa) 

c 95% 

t l 

c 70% 

t l 

Ff (kPa) 

(cl -o-3) u(t 
(kPa) 

R f 

100 

157.4 

149.53 

110.18 

0.0041 

0.0011 

266471.35 

172.87 

0.9104 

200 

294.4 

279.68 

206.08 

0.0093 

0.0048 

77223.283 

455.05 

0.6469 

300 

296.6 

281.77 

207.62 

0.0125 

0.0039 

138664.26 

336.26 

0.8820 


In the triaxial testing, soil is under 3 containment pressure levels. Each of these values correspond 
to parameter o 3eh . We determine the parameter E eiq . 


2.3.5. Determination of E e i q 

The formula to estimate Young’s modulus and its uncertainty involving the term 


Eeiq — 


ZiL ifjOi - 


(41) 


EitiW) 

To determine the range of elasticity, we will do calculations on the N°1 sample and apply it to all 
other samples. 
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For <j 3 — 100 kPa 


(°1 — a l)eh — 


2 • [0.23 • cos( 27.1) + 0.54 • sm(27.1)] 
1 1 — co.s(27. 1 ) | 


x 100 kPa 


Oi - Os) e h = 161,22 kPa 
The range of elasticity is given by the formula: 


^eh 


(°1 — °3 )e/i 



(°~1 q ~3)e/i 

(°1 ~ a 3)ult 


161,22 

266471,352 (l _ 22HL) 

V 229,523/ 


0,0020 


After having determined the three values of Young’s modulus at different confinement pressures, 
the value retained for a sample will be the average of the three. Calculations performed on our No. 1 
sample are presented in Table 4: 


Table 4 Estimation of Young’s modulus values for the B1N1 A WARAK site (Sample N°l) 


a 3 (kPa) 

£ eh 

Eelqi (KPa) 




A E elq \(kPa) 

E e iq(MPa ) 

Eelq 

100 

0.0020 

77223.2837 

0.0329 

879.1277 

171.45 

200 

0.0028 

171847.318 

0.1138 

300 

0.00125 

266471.352 

0.0633 


2.3.6. Determination of parametersE secl/2 and E sec 1/3 

Both parameters were calculated using formulas presented in paragraph 2.2.3. Table 5 presents the 
results of calculation performed on the sample N°1 (BINI A WARAK). 


Table 5 Young’s moduli values calculated on the B1N1 A WARAK samples (sample N°l) 


E° t ( Kpa) 

Rf 

Esec 1/2 

(Kpa) 

Esec 1/3 

(Kpa) 

E secl /2 

(Mpa) 

E secl/3 

(Mpa) 

266471.35 

0.9104 

145163.45 

185599.42 

91.63 

114.68 

77223.28 

0.6469 

52243.191 

60569.885 

138664.27 

0.8820 

77510.41 

97895.031 


3. RESULTS OF OUR APPROACH 


We present here four categories of results. The first concerns results of soil identification tests 
(Dmax> PI> Dso >f> S’ L,A) as well as the elastic -plastic parameters (Table 6); the second category of 
results represents the values of Young moduli calculated based on multiple correlations with soil 
identification parameters (Table 7). There will therefore be E LCPC calculated with the first correlation 
and the E SIKALI calculated with the second. The third category of results relates to the Young's moduli 
values deduced from triaxial test results based on the hyperbolic model of Duncan and Wong (Table 
8 ). 

We will have: 


fi) 0 : initial tangent modulus, 
fre^The equivalent modulus, 

E seci/ 3 et Eseci/2- The secant moduli. 
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The results of the fourth category allow us to estimate the correspondence between the two 
categories of Young’s moduli values, we have drawn the lines of linear regressions and calculated the 
correlation coefficients between each value of the second and third group. The figures 9 to 16 and the 
tables 9 to 10 summarize the expressions of straight linear regressions obtained and table 1 1 highlights 
the most relevant matches. 


Table 6 Identification of soil parameters 


Sites 

N° 

Name 

D„, ax (mm) 

Dso(mm) 

PI(%) 

G(%) 

S(%) 

L(%) 

A (%) 

C'(bars) 

(J)'(degree) 


1 

PG101 

20 

2.5 

29.3 

49.25 

29.28 

6.38 

14.11 

0.23 

27.1 


2 

PG102 

10 

4.5 

38.6 

49.56 

17.05 

15.94 

17.45 

0.15 

32 


3 

PG103a 

12.5 

0.4 

20 

41.49 

34.19 

8.57 

15.75 

0.49 

30.9 

* 

< 

4 

PG103b 

10 

2 

20 

38.61 

30.78 

14.55 

16.06 

0.53 

21.1 

K 

< 

5 

PG104a 

10 

0.7 

19.7 

33.81 

47.16 

7.98 

11.05 

0.08 

35.4 

£ 

6 

PG104b 

15 

9 

15.1 

40.365 

30.825 

11.585 

17.225 

0.13 

34.9 

< 

7 

PDlOla 

12.5 

2.6 

22.4 

55 

27.37 

6.35 

11.21 

0.17 

32.1 

s 

8 

PDlOlb 

10 

2.5 

30 

46.18 

27.95 

7.44 

11 

0.38 

27.9 


9 

PD 102 

10 

0.5 

20.5 

36.99 

31.2 

16.81 

15 

0.1 

30.5 


10 

PD 103a 

12.5 

7.5 

22 

60.79 

16.99 

13.21 

9.01 

0.03 

37.3 


11 

PD 103b 

15 

4 

15.1 

40.37 

30.83 

11.59 

15.23 

0.45 

28 


12 

PG101 

12.5 

0.005 

13 

35 

25 

20 

20 

0.36 

28 


13 

PG102 

10 

0.01 

14 

45 

21 

19 

15 

0.36 

29 


14 

PG103a 

12.5 

0.5 

15 

45 

22 

13 

20 

0.31 

39 


15 

PG103b 

15 

0.045 

29 

35 

33 

10 

22 

0.4 

31 


16 

PG104a 

12.5 

0.012 

21 

46 

26 

10 

18 

0.28 

37 

<& 

17 

PG104b 

12.5 

0.0675 

20 

32 

28 

21 

19 

0.75 

30 

o 

z 

18 

PDlOla 

5 

0.29 

25 

35 

32.5 

21 

11.5 

0.32 

34 

< 

Oh 

19 

PDlOlb 

15 

0.29 

25 

31 

29 

25 

15 

0.09 

21 

5 

20 

PD 102 

12.5 

0.014 

19 

33 

20 

25 

20 

0.28 

33 

o 

-J 

21 

PD 103a 

10 

0.05 

25 

36 

30 

19 

15 

0.38 

33 


22 

PD 103b 

15 

1.2857143 

34 

35 

25 

3 

37 

0.5 

29 


23 

PD 104a 

12.5 

0.017 

14 

36 

27 

19 

18 

0.16 

32 


24 

PD 104b 

12.5 

0.022 

15 

34 

32 

14 

20 

0.45 

30 


25 

PEG204b 

12.5 

4.8 

20 

27 

31 

12 

30 

0.2 

35 


26 

PED303a 

15 

0.01 

14 

47 

25 

13 

15 

0.75 

30 
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Table 7Young’s moduli values calculated on the basis of the correlations 


Sites 

N° 

Calculated Young’s moduli 

Elcpc (MPa) 

Esikali (MPa) 


1 

153.26 

81.34 


2 

63.05 

38.21 

X 

3 

94.15 

52.36 

i 

4 

84.49 

49.85 

hH 

< 

5 

90.95 

50.32 

£ 

6 

149.68 

87.61 

< 

7 

118.71 

66.69 

Z 

8 

81.711 

46.78 

3 

9 

76.64 

43.46 


10 

129.91 

73.04 


11 

145.56 

82.73 


12 

62.57 

33.21 


13 

60.17 

32.52 


14 

88.62 

51.92 


15 

68.70 

36.53 

o t 

16 

64.22 

34.08 

< 

T N 

17 

68.45 

37.21 

U 

z 

18 

39.20 

23.36 

< 

Clh 

19 

93.14 

48.68 

% 

20 

53.94 

29.05 

o 

21 

58.74 

31.60 

— 1 

22 

54.82 

34.57 


23 

72.66 

38.80 


24 

70.64 

38.28 


25 

74.22 

46.84 


26 

93.37 

47.48 


Table 8 Values of the Young’smoduli inferred from the results of the triaxial test 


Sites 

N° 

Experimental moduli 

£?( MPa) 

E efq (MPa) 

E sec i/ 2 (MPa) 

Eseci/ 3 (MPa) 


1 

160.79 

173.48 

91.63 

114.68 


2 

49.67 

82.62 

30.58 

36.95 


3 

93.67 

99.25 

53.28 

66.752 

t 

4 

66.34 

77 

46.72 

53.26 

< 

5 

73.06 

80.88 

46.97 

55.67 

£ 

6 

159.19 

169.58 

89.41 

112.67 

< 

7 

113.16 

115.13 

65.92 

81.66 

Z 

8 

90.38 

97.58 

55.39 

67.06 

s 

9 

73.16 

68.64 

46.93 

55.67 


10 

137.96 

147.26 

77.75 

97.82 


11 

153.33 

168.43 

83.35 

106.68 

D£ 

§ < 

12 

59.52 

56.042 

32.94 

40.64 

13 

59.06 

62.19 

39.21 

44.43 

|g 

14 

77.39 

78.51 

16.39 

19.29 

Oh 

15 

224.88 

226.93 

118.56 

154.00 

16 

79.13 

61.81 

24.93 

27.67 
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17 

62.62 

63.14 

24.93 

27.67 


18 

43.76 

40.95 

28.24 

32.47 


19 

97.67 

89.21 

67.22 

74.55 


20 

99.40 

100.07 

68.36 

78.71 


21 

66.99 

58.33 

31.88 

36.14 


22 

324.06 

285.15 

153.35 

197.29 


23 

89.16 

91.2 

58.24 

68.55 


24 

138.5 

139.17 

75.08 

96.22 


25 

64.06 

67.97 

37.01 

47.14 


26 

67.973 

120.83 

42.89 

51.25 


Analysis of the scattered values between experimental and calculated moduli highlights the 
presence of two outliers that were eliminated following the recommendations of the analytical 
statistics. We then traced the curves trends on the remaining scattered values. These trends curves are 
based on the least squares method (linear model). 


E° = f (ELCPC) 

180 


y * 1,0303 m 
R* =0.7895 



O ♦ 

O 20 40 60 80 100 120 140 160 180 


eicpc (MPa') 


Figure 9 Trend curve between and E LC pc 
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Esecl/2=f(ELCPC) 
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Figure 11 Trend curve between E secl / 2 and E LC pc 

Esecl/3=f(ELCPC) 


y=0,713x 
R 2 = 0,7057 


• t 


20 40 60 80 100 120 140 160 

ELCPC (MPa) 


180 

160 

40 

20 

00 

80 

60 

40 

20 


Figure 12 Trend curve between E secl / 3 and E LC pc 

Ef = f (ESI kali ) 
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Table 9 Summary of the regression lines depending on E LC pc 


Young modulus 

Linear regressions 

R-Square values 

£V° = f (ELCPC) 

y = 1.0303x 

R 2 =0.7895 

E elq = f (ELCPC) 

y = 1.1021x 

R 2 =0.8334 

Esecl/2 = f (ELCPC) 

y = 0.5838x 

R 2 =0.6866 

Esecl/2 = f (ELCPC) 

y = 0.713x 

R 2 =0.7057 


Table 10 Summary of the regression lines depending on Esikali 


Young modulus 

Linear regressions 

R-Square values 

£)° = f (ESIKALI) 

y = 1.8377X 

R 2 =0.7613 

E elq = f (ESIKALI) 

y = 1.9656x 

R 2 =0.8057 

Eseci/2 = f (ESIKALI) 

y = 1.0403x 

R 2 =0.6531 

Eseci/2 = f (ESIKALI) 

y = 1.2713x 

R 2 =0.6799 


Table 11 Selected regression lines for our local data 


Young modulus 

Linear regressions 

R-Square values 

Ef = f (ELCPC) 

y = 1.0303x 

R 2 =0.7895 

E elq = f (ELCPC) 

y = 1.1021x 

R 2 =0.8334 

Eseci/2 = f (ESIKALI) 

y = 1.0403x 

R 2 =0.6531 


4. CONCLUSION 

For pavement design some authors use the initial tangent modulus, others the secant moduli at 
predetermined points. The theory of soil mechanics defines several Young’s moduli and empirical 
calculation formulas of these parameters do not often specify the nature of the modulus. Then it is 
important to specify the nature of the Young’s modulus that our correlations describe through the 
evaluation of the regression lines between calculated and measured values (Table 9 and 10). This 
evaluation between estimated and calculated Young modulus has been obtained with our collected 
samples and the regression equations have provided significative results with a high R-Square value 
(Table 1 1). We saw that the LCPC modulus is closer to the initial tangent modulus and the equivalent 
modulus. The Sikali modulus approaches the secant modulus at 50% of the breaking stress. However, 
considering the number the relatively small samples, these results need to be confirmed or disproved 
by expanding further bank data necessary for the pursuit of this study. 
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